Separability Properties of Three-mode Gaussian States 
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We derive a necessary and sufficient condition for the separability of tripartite three mode Gaus- 
sian states, that is easy to check for any such state. We give a classification of the separability 
properties of those systems and show how to determine for any state to which class it belongs. We 
show that there exist genuinely tripartite bound entangled states and point out how to construct 
and prepare such states. 
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I. INTRODUCTION 

Entanglement of composite quantum systems is central 
to both the pecuharities and promises of quantum infor- 
mation. Consequently, the study of entanglement of bl- 
and multipartite systems has been the focus of research in 
quantum information theory. While pure state entangle- 
ment is fairly well understood, there are still many open 
questions related to the general case of mixed states. The 
furthest progress has been made in the study of systems 
of two qubits: it has been shown that a state of two 
qubits is separable if and only if its partial transpose is 
positive (PPT-property) [|| and a closed expression for 
the entanglement of formation was derived Moreover, 
it was shown that all entangled states of two qubits 
can be distilled into maximally entangled pure states by 
local operations. This property of distillability is of great 
practical importance, since only the distillable states are 
useful for certain applications such as long-distance quan- 
tum communication, quantum teleportation or cryptog- 
raphy g. 

In higher dimensions much less is known: the PPT- 
property is no longer sufficient for separability as proved 
by the existence of PPT entangled states in (D^ (8) (D* 
systems These states were later shown to be bound 
entangled even if two parties (Alice and Bob) share 
an arbitrarily large supply of such states, they cannot 
transform ("distill") it into even a single pure entangled 
state by local quantum operations and classical commu- 
nication. Meanwhile, a number of additional necessary 
or sufficient conditions for inseparability have been found 
for finite dimensional bipartite systems, which use prop- 
erties of the range and kernel of the density matrix p and 
its partial transpose p^^ to establish separability {M and 
references therein). 

When going from two to more parties, current knowl- 
edge is even more limited. Pure multipartite entangle- 
ment was first considered in ||^. A classification of TV- 
partite mixed states according to their separability prop- 
erties has been given . But even for three qubits there 
is currently no general way to decide to which of these 
classes a given state belongs [|lO|. Results on bound 



Recently increasing attention was paid to infinite di- 
mensional systems, the so-called continuous quantum 
variables (CV), in particular since the experimental re- 
alization of CV quantum teleportation Quan- 
tum information with CV in general is mainly con- 
cerned with the family of Gaussian states, since these 
comprise essentially all the experimentally realizable CV 
states. A practical advantage of CV systems is the 
relative ease with which entangled states can be gen- 
erated in the lab First results on separabil- 
ity and distillability of Gaussian states were reported in 
|l6|, ^ |20[ 1^, H. One finds striking similari- 



entanglement [|ll| and entanglement distillation |12 for 
multi-party systems have been obtained. 



^ ^ |20l |21|, |2S 

ties between the situations of two qubits and two one- 
mode CV systems in a Gaussian state: PPT is necessary 
and sufficient for separability jl7[ |l^ , and all inseparable 
states are distillable Generalizing the methods re- 

viewed in 1^] it was shown that for more than two modes 
at either side PPT entangled states exist In a 
computable measure of entanglement for bipartite Gaus- 
sian states was derived. 



The study of CV multipartite entanglement was initi- 
ated in [in 1^ , where a scheme was suggested to create 
pure CV A'^-party entanglement using squeezed light and 
A^— 1 beamsplitters. In fact, this discussion indicates that 
tripartite entanglement has already been created (though 
not investigated or detected) in the CV quantum telepor- 
tation experiment |14|. 



In this paper we provide a complete classification of 
tri-mode entanglement (according to the scheme |9j ) and 
obtain - in contrast to the finite dimensional case - a sim- 
ple, directly computable criterion that allows to deter- 
mine to which class a given state belongs. We show that 
none of these classes are empty and in particular provide 
examples of genuine tripartite bound entangled states, 
i.e. states of three modes A, B, and C that are separable 
whenever two parties are grouped together but cannot be 
written as a mixture of tripartite product states. 

Before we can derive our results we need to introduce 
some notation and collect a number of useful facts about 
our main object of study: Gaussian states. 



2 



II. GAUSSIAN STATES 

In quantum optics and in other scenarios described by 
continuous quantum variables, not all states on the infi- 
nite dimensional Hilbert space are equally accessible in 
current experiments. In fact, the set of Gaussian states 
comprises essentially all genuinely CV states that can 
currently be prepared in the lab. This, and the mathe- 
matical simplicity of these states are the reasons why CV 
quantum information has so far considered almost exclu- 
sively Gaussian states, as will the present paper. This 
section summarizes results on Gaussian states that we 
need in the following and introduces some notation. 

We consider systems composed of n distinguishable in- 
finite dimensional subsystems, each with Hilbert space 
Ti, = L^(]R). These could be implemented quantum op- 
tically by different modes of the electromagnetic field, 
hence each of these subsystems will be referred to as a 
"mode". To each mode belong the two canonical ob- 
servables X^, P^, fc = 1, . . . , n with commutation relation 
[Xk.Pk] = i- Defining Rk = Xk,Rn+k = Pk these rela- 
tion are summarized as [Rk,Ri\ = —iJki, using the anti- 
symmetric 2n X 2n matrix 



J ^ 



-1 

1 



(1) 



which plays an important role in the following calcula- 
tions II . 

For such systems, it is convenient to describe the state 
p by its characteristic function 

X{x) = ii[pD{x)]. (2) 

Here x — {q,p), q,p G R" is a real vector, and 

D{x) = e-''S'=(«'=-^'=+P''-f''=). (3) 

The characteristic function contains all the information 
about the state of the system, that is, one can construct 
p knowing x. Gaussian states are exactly those for which 
X is a Gaussian function of the phase space coordinates 
X M, 



(4) 



where 7 is a real, symmetric, strictly positive matrix, the 
correlation matrix (CM), and d G IR^" is a real vector, 
the displacement. Note that both 7 and d are directly 
measurable quantities, their elements 7^; and dk are re- 
lated to the expectation values and variances of the oper- 
ators Rk- A Gaussian state is completely determined by 
7 and d. Note that the displacement of a (known) state 
can always be adjusted to d = by a sequence of unitaries 
applied to individual modes. This implies that d is irrel- 
evant for the study of nonlocal properties. Therefore we 
will occasionally say, e.g., that "a CM is separable" when 
the Gaussian state with this CM is separable. Also, from 
now on in this paper "state" will always mean "Gaussian 
state" (unless stated otherwise). 



Not all real, symmetric, positive matrices 7 correspond 
to the CM of a physical state. There are a number of 
equivalent ways to characterize physical CMs, which will 
all be useful in the following. We collect them in 

Lemma 1 (Correlation Matrices) 

For a real, symmetric 2n x 2n matrix 7 > the following 
statements are equivalent: 



7 is the CM of a physical state, 



Jl'^J > 0, 



7-iJ> 0, 



(5a) 



(5b) 



(5c) 



-f = S^{D®D)S, (5d) 
for S symplectic Ml] and D > 1 diagonal \2& 



Proof 



^) ^ (^ see |6[ 
1) <^ see [pi Prop. 4.22]. 



( 



^ d5 



see 
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A CM corresponds to a pure state if and only if (iff) 
D = 1, i.e. iff det7 = 1 (e.g. p6|). It is easy to see from 
( ^ ) that for pure states Ineq. (|5b|) becomes an equality 
and dim[ker(7 — iJ)] — n. It is clear from Eq. ( |5d[ ) that 
for every CM 7 there exists a pure CM 70 such that 
7o ^ 7- This will allow us to restrict many proofs to 
pure CMs only. Note that for a pure 2n x 2n CM 7 it 
holds that tr7 > 2n. 

A very important transformation for the study of en- 
tanglement is partial transposition Transposition is 
an example of a positive but not completely positive map 
and therefore may reveal entanglement when applied to 
part of an entangled system. On phase space transposi- 
tion corresponds to the transformation that changes the 
sign of all the p coordinates {q,p) A((7,p) = {q,—p) 
p8[ and leaves the q's unchanged. For 7 and d this means 
(tVo?) i-^ (A7A, Ad). Using this, the NPT-criterion for in- 
separability 1^ translates very nicely to Gaussian states. 
Consider a bipartite system consisting of m modes on 
Alice's side and n modes on Bob's (m x n-system in the 
following) . Let 7 be the CM of a Gaussian m x n-state 
and denote by A^ = A © 1 the partial transposition in 
A's system only. Then we have the following criterion for 
inseparability: 

Theorem 1 (NPT criterion) 

Let 7 he the CM of a Ixn system, then 7 corresponds to 
a inseparable state if and only if AaJ-^a is not a physical 
CM, i. e. if and only if 



We say that 7 "is NPT" if ^ holds. 



(6) 
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Proof: See flj] for TV = 1 and [gO| for the general 
case. 

Occasionally it is convenient to apply the orthogonal 
operation to the right hand side of Ineq. and write 

Ja = I^aJ^a- 

For states of at least two modes at both sides Condi- 
tion (|^) is still sufhcient for inseparability, but no longer 
necessary as shown by Werner and Wolf, who have con- 
sidered a family of 2 x 2 entangled states with positive 
partial transpose pO| ]. In the same paper, it was shown 
that 



Theorem 2 (Separability of Gaussian States) 
A state with CM"f is separable iff there exist CMs ^a,Jb 
such that 



7 > 7^ © 7_B. 



(7) 



It is observed in [g^l that if Ineq. (0) can be fulfilled, 
then the state with CM 7 can be obtained by local op- 
erations and classical communication from the product 
state with CM 7p — 7a ffi 7b, namely by mixing the 
states {'jpjd) with the d's distributed according to the 
Gaussian distribution oc cx'p\—(f"{j — 'yp)~^d\. 

Note that while Theorem g gives a necessary and suffi- 
cient condition for separability, it is not a practical crite- 
rion, since to use it, we have to prove the existence or non- 
existence of CMs 7a, 7b. Instead, a criterion would allow 
to directly calculate from 7 whether the corresponding 
state is separable or not. Theorem |^ and its extension 
to the 3-party situation are the starting point for the 
derivation of such a criterion for the case of three-mode 
three-party states in the following main section of this 
paper. 



III. TRI-MODE ENTANGLEMENT 

When systems that are composed of > 2 parties are 
considered, there are many "types" of entanglement due 
to the many ways in which the different subsystems may 
be entangled with each other. We will use the scheme in- 
troduced in 1^ , to classify three-mode tripartite Gaussian 
states. The important point is that from the extension of 
Theorem |^ we can derive a simple criterion that allows to 
determine which class a given state belongs to. This is in 
contrast to the situation for three qubits, where up until 
now no such criterion is known. In particular, we show 
that none of these classes are empty and we provide an 
example of a genuine tripartite bound entangled state, 
i.e. a state of three modes A, B, and C that is separa- 
ble whenever two parties are grouped together but can- 
not be written as a mixture of tripartite product states 
and therefore cannot be prepared by local operations and 
classical communication of three separate parties. 



A. Classification 

The scheme of [|j considers all possible ways to group 
the N parties into m < N subsets, which are then them- 
selves considered each as a single party. Now, it has to 
be determined whether the resulting m-party state can 
be written as a mixture of m-party product states. The 
complete record of the m-separability of all these states 
then characterizes the entanglement of the A^-party state. 

For tripartite systems, we need to consider four 
cases, namely the three bipartite cases in which AB, 
AC, or BC are grouped together, respectively, and the 
tripartite case in which all A, B, and C are separate. 
We formulate a simple extension to Theorem |^ to 
characterize mixtures of tripartite product states 

Theorem ^' (Three-party Separability) 
A Gaussian three-party state with GM 7 can be written 
as a mixture of tripartite product states iff there exist 
one-mode correlation matrices ^a^IBtIC such that 



7 - 7A © 7s © 7c > 0. 



(8) 



Such a state will be called fully separable. 

Proof: The proof is in complete analogy with that of 

Theorem |^ in |po| and is therefore omitted here. 

A state for which there are a one-mode CM 7^ and 
a two-mode CM "fBC such that 7 — 7a © 7bc > is 
called A—BC biseparable (and similarly for the two other 
bipartite groupings). In total, we have the following five 
different entanglement classes: 

Class 1 Fully inseparable states are those which are not 
separable for any grouping of the parties. 

Class 2 1-mode biseparable states are those which are 
separable if two of the parties are grouped together, 
but inseparable with respect to the other groupings. 

Class 3 2-mode biseparable states are separable with re- 
spect to two of the three bipartite splits but insep- 
arable with respect to the third. 

Class 4 3-mode biseparable states separable with respect 
to all three bipartite splits but cannot be written 
as a mixture of tripartite product states. 

Class 5 The fully separable states can be written as a 
mixture of tripartite product states. 

Examples for Class 1 (the GHZ-like states of p^). 
Class 2 (two-mode squeezed vacuum in the first two and 
the vacuum in the third mode), and Class 5 (vacuum 
state in all three modes) are readily given; we will pro- 
vide examples for Classes 3 and 4 in Subsection IV below. 



How can we determine to which Class a given state 
with CM 7 belongs? States belonging to Classes 1, 2, 
or 3 can be readily identified using the NPT-criterion 
(Theorem |^) . Denoting the partially transposed CM by 
7a; — Kx^P^x, X — A,B, C, we have the following equiva- 
lences: 
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Lemma 2 ( Classification ) 

lA i iJ,lB i iJ,lc ^iJ ^ Class 1 (9) 

{*)lA t iJ.lB t iJnc >iJ ^ Class 2 (10) 

{*)lA ^ iJ.lB > iJilc >iJ^ Class 3 (11) 

lA > iJ,lB > iJ,lc >iJ'^ Class 4 or 5, (12) 

where the (*) reminds us to consider all permutations of 
the indices A, B, and C. 

The proof follows directly from the definitions of the dif- 
ferent classes and Theorem ||. 

What is still missing is an easy way to distinguish be- 
tween Class 4 and Class 5. Thus to complete the classifi- 
cation we now provide a criterion to determine whether a 
CM 7 satisfying Ineqs. (12) is fully separable or 3- mode 
biseparable, that is we have to decide whether there exist 
one- mode CMs 7a, 7b, 7c such that (^) holds, in which 
case 7 is fully separable. In the next subsection we will 
describe a small set consisting of no more than nine CMs 
among which 7^ is necessarily found if the state is sepa- 
rable. 



B. Criterion for Full Separability 

This subsection contains the main result of the pa- 
per: a separability criterion for PPT 1x1x1 Gaussian 
states, i.e. states whose CM fulfills Ineqs. (^2|). We start 
from Theorem ^' and obtain in several steps a simple, 
directly computable necessary and sufficient condition. 
The reader mainly interested in this result may go di- 
rectly to Theorem ^, from where she will be guided to 
the necessary definitions and Lemmas. 

Since the separability condition in Theorem is for- 
mulated in terms of the positivity of certain matrices the 
following lemma will be very useful throughout the pa- 
per. We consider a self-adjoint {n + m) x (rt + m) matrix 
M that we write in block form as 



M = 



A C 
Ct B 



(13) 



where A^ B, C are n x n,m x m, and n x m matrices, 
respectively. 

Lemma 3 (Positivity of self-adjoint matrices) 
A self-adjoint matrix M as in with A > 0, B > is 
positive if and only if for all e > Q 

A-C^C^>0, (14) 

or, equivalently, if and only if 

kerSCkerC (15a) 

and 

A-C^C^>Q, (15b) 
B 

where B~^ is understood in the sense of a pseudoinverse 
(inversion on the range). 



Proof: The only difficulty in the proof arises if ker B ^ 
0. Therefore we consider the matrices Mg, where B in 
( |T3| ) is replaced by i?e — B + el (e > 0), which avoid this 
problem and which arc positive Ve > iff M > 0. In a 
second simplifying step we note that > Ve > iff 
M'^ = (1 e B7^'^)M{1 e 57^^^) > 0. 

Now direct calculation shows the claim: we can write a 



general / © 5 as / ( 



where is or- 



(S7'/'ct)/i + / 

thogonal to the range of (Be^^^^Ct). Then {f®g)^M'^{f® 

g) - f\A-CB-^C^)f+{f + h;)^CB-'C\f + h) + h\hu 
which is clearly positive, if (O) holds. With the choice 
h± — G and h — —f ii is seen that (|4|) is also necessary. 

That the second condition is equivalent is seen as fol- 
lows: If Ineq. ( p^ holds Ve > there cannot be vec- 
tor ^ G kerB and ^ ^ kerC since for such a ^ we have 

Cm C < for sufficiently small e > 0, 



3 + £l 



and if ( p^ h olds' then (|l|) converges to ( |l5b[ ). Con- 
versel y, if ( |15a| ) holds, then CB~^C^ is well-defined and 
Ineq. ( |15b| ) imphes it Ve > 0. ■ 
As mentioned above, in this section we exclusively con- 
sider three-mode CMs 7 that satisfy Ineqs. (p^). We 
write 7 in the form of Eq. ( p^ ) as 



A 



(16) 



where ^ is a 2 x 2 matrix, whereas i? is a 4 x 4 matrix. 
We observe that Ineqs. (|l2| ) impose some conditions on 
7 that will be useful later on: 



(17) 



Observation 1 Let 7 satisfy Ineqs. (jl^), then 

(TAiJ 
7 > I (TBiJ 
aciJ 

where cr^ £ {0,±l},Va; = A,B,C. 



Proof: Ineqs. < \12\ j say that j±iJ >0 and j±iJx > 
Vx. By adding these positive matrices all combinations 
of (Tj. can be obtained. ■ 
From this it follows: 



Observation 2 For a PPT CM 7 as in Eq. I\U 

ker(B + iJ), ker(B -I- iJ) C ker C, (18) 

where J = J ® (—J) is the partially transposed J for two 
modes. 

Proof: Cond. ([l^) on the kernels is an immedi- 
ate consequence of Lemma || applied to the matrices 
7 — ® I J ® (±« J), which are positive by Obs. |l|. ■ 
Then the matrices 



N 



A-C- 



1 



B-iJ 



N = A-C- -C^ 

B — iJ 



(19a) 
(19b) 



are well-defined and 
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Observation 3 It holds that both 

trN,tTN>0, (20) 

Proof: Cond. (|20| ) is true since, again by Lemma ^ and 
Obs. 0, both N and Nave positive and N±iJ, N±iJ > 0. 
This imphes that N, N cannot be zero, which is the only 
positive matrix with vanishing trace. Therefore trA^, trA'^ 
are strictly positive. ■ 
The remainder of this section leads in several steps to 
the separability criterion. First, we simplify the condition 
® by reducing it to a condition which involves only one 
one-mode CM 7^. 

Lemma 4 A PPT 3-mode CM 7 is fully separable if and 
only if there exists a one-mode CM 7^ such that 



N > 7A, 
N > 7A, 



(21a) 
(21b) 



where N, N were defined in Eqs. Without loss 

of generality we require to be a pure state CM, i.e. 
det7^ = 1. 

Proof: By Theorem H' full separability of 7 is equiv- 
alent to the existence of one-mode CMs 7^, 7b, 7c > iJ 
such that 7 — Ta ffi 7b ffi 7c > 0. Let stand for 7a,b.c- 
By Lemma this is equivalent to 37^; such that = 



C >^B® ic, Ve > 0, where A^ = A + el. 
But iff there exist such 7^, then (Lemma ^) th e ine quality 
also holds for e = and the kernels fulfill (15a). This 
is true iff the matrix X = Xq is a CM belonging to a 
separable state, i.e. (Theorem ||) iff X' > iJ,iJ. Using 
B > iJ, iJ [which holds since 7 fulfills Ineqs. (^2|)] we 
obtain that 7 is separable iff there exists 7_a > such 
that 



A -J A 



C 
Bi 



> 0, A: = 1,2, 



(22) 



wher e B'^ = B — iJ and B'2 — B — iJ. Since Condition 



(15a) holds, this is (Lemma ^ equivalent to Ineqs. (21). 
That we can always choose det7yi = 1 follows directly 

from Eq. ( [5d| ) and the remark after Lemma |l|. m 

While we can always find a ja fulfilling Ineq. (21b), 



since 7 belongs to a PPT state (and there exists a two- 
mode CM jBC > iJ such that 7^ ©tbc is smaller than 
7), it may well happen that Ineq. ( plaj ) cannot be satis- 
fied at all, or that it is impossible to have both Ineqs. ( |2l| ) 
fulfilled for one ja simultaneously. Note that due to In- 
eqs. ( p^, N and N as above are always positive. From 
Ineqs r^pi]) we observe that 

Observation 4 for the CM j of a separable state it is 
necessary to have 



tiN,tiN > 2, 
detA^,detiV > 0, 



(23a) 
(23b) 



where 7 as in Eq. ( |7^ j and N as in Eqs. (|7 



Proof: A self-adjoint 2x2 matrix is positive iff its trace 
and determinant are positive. Since the trace of the RHS 
of both Ineqs. ( |2l| ) is > 2 [remark after Lemma |^ the 
same is necessary for the LHS. Also, since det7^ — 1, 
which implies that ja has full rank, any matrix > ja 
must also have full rank ]3l|| and thus a strictly positive 
determinant. ■ 
For a self-adjoint positive 2x2 matrix 



R 



a b 
b* c 



(24) 



we show 



Lemma 5 There exists a CM ja l£ R o,nd only if there 
exist (y, z) £ IR^ such that 



trR > 2y/l + y^ + z^, (25a) 



det i? + 1 + ( ^ ) > tri? VTTy^T^, (25b) 



whe 



L= {a~c, 2Reb). 



(26) 



Proof: As noted in Lemma ^ we need only look for 7^ 
with detjA — 1- We parameterize 



lA 



X + y z 
z X ^ y 



(27) 



with real parameters x, y, z and x^ — 1 + + for 
purity. This is a CM iff 7a — > (Lemma |l|), that 
is iff ii^A = 2a; > (where we use that positivity of the 
a 2 X 2 matrix is equivalent to the positivity of its trace 
and determinant and det(7A — i J) = by construction). 
By the same argument, i? — 7^ > leads to the two 
conditions (p5[). ■ 

The inequalities (|2^) have a simple geometrical inter- 
pretation that will be useful for the proof of the promised 
criterion: Ineq. (25a) restricts (y, 2;) to a circular disk 
C of radiu s (tri?)-^/4 — 1 around the origin, while 
Ineq. ( |25bD describes a (potentially degenerate) ellipse 
£ (see Fig. |), whose elements are calculated below, and 
the existence of a joint solution to Ineqs. ( p5| ) is therefore 
equivalent to a nonempty intersection of C and £. 

Applying this now to the matrices (|l^) we find that 
in order to simultaneously satisfy both conditions in 
Lemma ^, the intersection between the two ellipses £,£ 
and the smaller of the two concentric circles C',C' (which 
we denote in the following by C) must be nonempty. This 
condition leads to three inequalities in the coefficients 
of the matrices TV, A^ which can be satisfied simultane- 
ously if and only if the PPT trimode state is separable. 
Thus we can reformulate the condition for separability 
(Lemma |^) as follows 

Lemma 6 (Reformulated Separability Condition) 

A three-mode state with CM 7 satisfying Ineqs. ([Tlj is 
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fully separable if and only if there exists a point {y, z) S 
fulfilling the following inequalities: 



min{triV, triV} > 2y/\ + ipT^, (28a) 
AetN + l + L^(^^ > tiNy/l + y^ + 7J, (28b) 

det N + I + L'^ (^^ > trN^/lTy^^^. (28c) 

Proof: According to Lemma ^ 7 belongs to a separable 
state iff we can find smaller than N and smaller than 
N. According to Lemma ^ we can find such a iff we 
can find (y, z) such that Incqs. ( p5|) are satisfied for both 
N and N. u 

In the following paragraphs wc have a closer look at 
the sets 5,5, and C. The goal of this discussion is to 
identify a few special points - directly computable from 
7 - among which a solution to Ineqs. (|2^) will be found 
iff the state under consideration is separable. This will 
then lead to the final practical form of the separability 
criterion which is stat ed a t the end of this section. 

By squaring Ineq. ( ^8b| ) we obtain 



K 



< m. 



(29) 



where fi = (detiV + l)/ki, 
and the matrix K is [M 



[(det7V + l)2-fci], 



K = kiPL + k2PL^, 

with the orthogonal projectors Pl, Pl^ on L, and 

ki = A [detiV+ (Im6)2] , 
k2 - {tvNf. 

Due to Ineqs. ( p^ ) ki and k2 are strictly positive, ^, m are 
well-defined and K \s & positive matrix of rank 2. Let us 
now distinguish the cases m < and m > 0. For m < 
Ineq. (|2^) can never be fulfilled since X is a positive 
matrix. In the case m > 0, Ineq. (|2|) describes an ellipse 
£ which is centered at rrie — fiL with major axis L and 
minor axis L-^ of lengths m/ki > ^Jrnjhi^ respectively. 
From Ineq. (28c) we obtain the same equations for the 
tilded quantities derived from N . 

The final argument for the derivation of the separabil- 
ity criterion is as follows. By Lemma ^ the state is sepa- 
rable if and only if the three sets described by Ineqs. ( |2^ ) 
have a common intersection, i.e. iff/ = £n5nC7^0. 
The border of / is contained in the union of the borders 
of the ellipses and circle: dl C 9£U9£U9C Now we can 
distinguish two cases, both of which allow to calculate 
a definite solution to the Ineqs. (p8| ) if the state is sep- 
arable: Either dl has nonempty intersections with the 
borders of two of the sets £,£,C 01 dl coincides with the 
border of one of the three. In the latter case this whole 
set is contained in /. In the former case, at least one of 
the points at which the borders intersect must be in / 



and thus a solution. If no solution is found this way the 
state is inseparable. This argument is made more precise 
in the final theorem. Formulas for the nine candidate 
solutions - the centers mc, me,mg and the intersections 
points if-, ifg, if- - are given in Appendix 

Theorem 3 (Criterion for full separability) 
A three-mode state corresponding to the CM 7 satisfying 
Ineq. (|7Jj is fully separable if and only if Ineq. ^3l\ ) holds 
and there exists a point ^^oi 



s: 'ceJ 



(30) 



fulfilling the Ineqs. ([^ 



Proof: We already saw (Obs. §) that det N, det N > 
are necessary for separability. If this holds, the quanti- 
ties used in (|2|, ^o|) and in their derivation are all well- 
defined. 

According to Lemma ^ 7 is fully separable iff there ex- 
ists a point {y, z)"^ such that the Ineqs. ( p8| ) are fulfilled. 
Therefore, if one of the points (|3^) satisfies Ineqs. ( |2^ ) 
then it determines a 7^ fulfilling Ineqs. ( |2l| ) thus prov- 
ing that the state is separable. To complete the proof, we 
show that if the state is separable, then we find a solution 
to Ineqs. ( p8| ) among the points (|30|). 

As pointed out before, the condition that Ineqs. (|2^ ) 
can simultaneously be satisfied has the geometrical inter- 
pretation that the circle C and the two ellipses £,£ have 
a nonempty intersection, i.e. I = £ CiS CiC ^ ^. 

Thus it remains to prove that if / is nonempty then 
one of the nine points in (|3^) lies in /. But if / ^ there 
are only the following two possibilities: since all the sets 
considered are convex and closed, either the border of 
/ coincides with that of one of the sets C,£,£ (which 
means that one of these sets, call it 5, is contained in 
both others) or at least two of the borders dC,d£,d£ 
contribute to dl, in which case the points at which these 
two intersect belong to dl and thus to /. 

In the former case, the center of 5 is a solution and 



given by one of the Eqs. (Al); in the latter, one can find 
a solution among the intersections of the borders of the 
sets £,£,C. That these are given by the i^ is shown in 
Appendix ■ 
If a CM 7 belongs to a separable state according to the 
above theorem then the point ^soi provides us with a pure 
one-mode CM ja such that N, N > ja- By construction 
7' = B - C{A - ja)~^C^ is a separable 2 x 2 CM and 
by repeating a similar procedure as above with 7' we 
can calculate a pure product-state decomposition of the 
original state with CM 7. 



IV. EXAMPLES OF BOUND ENTANGLED 
STATES 

In this section we construct states belonging to Classes 
3 and 4. Our construction makes use of ideas that were 
first applied in finite dimensional quantum systems to 
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find PPT entangled states (PPTES) ||l and then gener- 
alized in [^2[ to construct so-called edge states, i.e. states 
on the border of the convex set of states with positive 
partial transpose. Similarly, one can define "edge CMs" 
as those that lie on the border of the convex set of PPT 



CMs (they are called "minimal PPT CMs" in fgO)). 

This section is divided into three subsections. In the 
first one we define "edge CMs" and characterize them. In 
the second and third subsections we present two different 
families of CMs which contain edge CMs. We also show 
that within those families we have CMs belonging to all 
classes. 



A. Edge CMs 

In the following we will consider CMs 7 corresponding 
to PPT states, i.e. fulfilhng 



-t-iJx> 0, for aU x = 0, ^, C, 



(31) 



where Jq = J. 

Definition 1 (Edge Correlation Matrices) 

A CM 7 is an edge CM if it corresponds to a non- 
separable state, fulfills ^Ij), and 7' = 7 — P does not 
fulfill d^jf ) for all real operators P with 7^ P > 0. 

Note that a state with an edge CM automatically be- 
longs to class 4 (i.e. edge CMs correspond to 3-mode 
biseparable states). In order to fully characterize them, 
we will need the following definition. Let us consider the 
complex vector space V Q oi dimension d spanned 
by the vectors belonging to the kernels of all j — iJx 
{x — 0, A, B,C). We will define K{'y) as a real vector 
space which is spanned by the real parts and imaginary 
parts of all the vectors belonging to V. More specifically, 
let us denote by B = {/^ + «// a basis of V, such 
that and Jj are real. We define 



imposing that there is no real vector in the intersection of 
the ranges of the matrices ^ — iJx- This is again equiva- 
lent to saying that there is no real vector orthogonal to all 
the ker(7 — iJx), which in turn is equivalent to K = R^, 
since that vector should be orthogonal to all the real and 
imaginary parts of the vectors spanned by those kernels. 
Now, if 7 corresponds to an entangled state it is clear 
that 7 7^ 7yi ffi 7_B ffi 7c- Conversely, if 7 7^ 7^ ® 7b ® 7c 
was separable, then there must exist some real positive 
P such that 7 — P = 7A ® 7b ® 7c is separable, and 
therefore fulfills (|3^), which is not possible. ■ 
Note that this Theorem generalizes easily to the cases 
of more than three parties and more than one mode at 
each site. 

In the construction of the following two examples of tri- 
partite bound entangled states we are going to use this 
theorem. The idea is to take a CM 70 of a pure entan- 
gled state [which, of course, does not fulfill (Pl|) ] and add 
real positive matrices until the conditions (|3l|) as well as 
K = R*^ are fulfilled. If the resulting CM is not of the 
form 7a ® 7b ® 7c then Theorem 4 implies that it is an 
edge CM. In fact, we can add more real positive matrices 
keeping the state entangled [and fulfilling (|3^)]. In order 
to see how much we can add, we can use the criterion 
derived in the previous section. 

This method of constructing CMs belonging to Class 4 
also indicates how the corresponding states may be pre- 
pared experimentally. Adding a positive matrix P to the 
CM 70 corresponds to the following preparation process: 
start with an ensemble of states with CM 70, displace 
them randomly by d according to the Gaussian proba- 
bility distribution with covariance matrix given by the 
inverse of P. This is a local operation (that potentially 
needs to be supplemented by classical communication) 
on each individual mode. The state produced by this 
randomization has CM 7 + P [M . 



B. Example 1 



X(7) 



:Afe/A + M/c//,Afe,MfeeR CR6. (32) 



that is, the real span of the vectors and fj. Note 
that this definition does not depend on the chosen basis 
B [As it is pointed out in Appendix ^ ^(7) coincides 
with the real vector space spanned by all the vectors in 
the kernels of 7 -I- Jxl~^Jx\- We then have the following 

Tiieorem 4 ( Characterization 0/ 1 x 1 x 1 edge CMs ) 
A CM 7 fulfilling (^Jp is an edge CM if and only if there 
exist no CMs ^a,1b,1c such that 7 = 7a® 7b® 7c aiT'd 
K = R6. 

Proof: We will use the fact that, given two posi- 
tive matrices A, P 7^ 0, there exists some e > such that 
A - eP > iff ran(P)C ran(A). According to the Def. |l| 
we cannot subtract any real positive matrix from 7 with- 
out violating the conditions (^TJ). This is equivalent to 



In the first example we start out with an entangled 
state between the two parties Alice and Bob and the vac- 
uum state in Charlie and add two projectors to the cor- 
responding CM. More specifically, we consider the CMs 
of the form 7ai.a2 = 7 + aiPi + a2-P2, where 



7 = lAB ® Ic, 



and 



7AB 



a c 

a -c 

c a 

-c a 



(33) 



(34) 



with a = Vl + and c can take any value different from 
zero. Here, Pi — pip[ and P2 — pipj , where pi — 
(0, 1, 0, 1, 1, 2)^ and P2 = (1, 0, -1, 0, 0, 1)^. 

In order to explain why the CM 7ai,a2 achieves our 
purposes, let us first consider the two-mode case in which 



the correlation matrix is jab- We denote now hy p = 
Pi + iP2 [where pi = (0, 1, 0, 1)^ and = (1, 0, -1, 0)^] 
the eigenvector corresponding to the negative eigenvalue 
of 7AS — iJa H^- Since (— iXi)* = —iJb we have that 
the eigenvector corresponding to the negative eigenvalue 
of JAB — iJB is p* = pi — ip2- By adding a sufficiently 
large multiple of the projectors onto those vectors, we 
obtain a CM whose partial transposes are positive. Note 
that in this case (just two modes) this would already 
make the state separable. 

In the case of three modes with a correlation matrix 7 
the same argumentation applies, namely that by adding 
some projectors we can make the partial transposes with 
respect to A and B positive. However, we have to in- 
volve C and thereby smear out the initial entanglement 
between A and B among all three parties. This is ex- 
actly what is achieved by adding the projectors Pi and 
P2. If we choose now, for instance, c = 0.3, oi = 1, and 
02 ~ 0.5531095, then one can show that the set if (701,02) 
defined as in Eq. (|3^) spans IR^. As mentioned at the end 
of the previous subsection, since the resulting CM is not 
of the form 7^ © 7s © 7c it corresponds to an edge CM. 

In Fig. |l| we illustrate to which class 701,02 belongs as 
a function of the parameters ai.2. In order to determine 
this, we have used the criterion derived in the previous 
section. It is worth noting that 701,02 never becomes 
separable. This follows from Theorem |^ and the fact 
that both TO = m = for all values of 01,2, as can be 
easily verified. This implies that the two ellipses [c.f. 
Ineq. (|29|)] are just two points [which coincide with the 
centers given in Eq. (Al)]. Thus, the only possibility that 
the circle and the two ellipses intersect is that the centers 
of the ellipses are the same and lie inside the circle. It is 
easy to show that for all values of ai and 02 the centers 
of the two ellipses are never the same. Thus the state 
corresponding to the CM 701,02 is never separable and 
is a PPTES for all values of ai , 02 for which the partial 
transposes are positive. 



C. Example 2 

Here we present a family of states which belong either 
to Class 1, 4, or 5. The states of this family are obtained 
from a pure GHZ-like state [Q by adding a multiple of 
the identity, i.e.. 



7a = 7 + al, 



(35) 



where 
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(36) 
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FIG. 1: The Entanglement Classes of 7ai,a2 



with a > 1 and 



6= i(5a- V9a2-8), 
c = i(a- \/9a2 - 8). 



(37) 
(38) 



For the following discussion, we pick a — 1.2. It is 
clear that for a = the state is fully inseparable, i.e., 
it belongs to Class 1, whereas for a > 1 the state will 
be fully separable (Class 5). We will show now that for 
cto < 01 < at, where ao « 0.29756 and ai ~ 0.31355 the 
state is biseparable and belongs therefore to Class 4. 

The CM 7q, is symmetric with respect to permutations 
between the parties, and therefore the negative eigenval- 
ues of the matrices 7 — iJx,x = A,B,C are the same. 
We denote its absolute value by ao ~ 0.29756. It is easy 
to determine the real and imaginary part of the corre- 
sponding eigenvectors. One finds that all those vectors 
are linearly independent. If we add now aol to 7 then 
all those vectors belong to K{jao) which immediately 
implies that ^^(7^0) = R^. Since 7^ 7A ©7b ® 7c we 
have that it is an edge CM. 

Let us now use Theorem 3 in order to determine ai. 
First of all, we show, independently of the discussion 
above that belongs to Class 4. In particular, we find 
that m = TO = [cf. Eq. (p9|)], which implies that there 
exists a solution to the Ineqs. ( p8| ) only if the centers of 
the two ellipses are the same and lie within the circle. 
Here one can also show that the two centers are not the 
same and so the state corresponding to the CM •jag is 
a PPTES. Let us determine the values of a for which it 
is still the case that there exists no intersection of the 
two ellipses and the circle given by the Ineqs. (p8[). It 
is easy to show that if a > ao then triV < triV, which 
implies that the circle that has to be considered has ra- 
dius Tc = y/ (trAf)2/4 — 1. One can also easily verify that 
the two ellipses never intersect the border of the circle. 
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FIG. 2: (a) The circle and the two eUipses do not have a 
joint intersection, therefore the state corresponding to 7a is a 
PPTES, (b) the circle and the two ellipses have a joint inter- 
section, therefore the state corresponding to 7c, is separable. 



which simplifies the problem. The ellipses must always 
lie inside the circle (since if they where outside it would 
never be possible to obtain a separable state even for 
a > 1). Thus, the problem reduces to check at which 
point the ellipses intersect each other. This occurs when 
a = ai w 0.31355. Thus the CM 7^, where < a < a\ 
corresponds to a PPTES, whereas for a > ai, the cor- 
responding state is fully separable. In Fig. ^ we have 
plotted the circle and the two ellipses, which are almost 
circles in this case, for (a) a < oi\ and (b) a> a\. 



V. CONCLUSIONS 

We have discussed nonlocal properties of Gaussian 
states of three tripartite modes. We have distinguished 
five classes with different separability properties and 
given a simple necessary and sufficient criterion that al- 
lows to determine which of these classes a given Gaus- 
sian state belongs to. The first three classes contain only 
NPT states and positivity of a state under the three par- 
tial transpositions suffices to determine to which of those 
it belongs. The separability criterion, which allows to 
distinguish PPT entangled states from separable states 
is the main result of this paper. For the case of three 
qubits such a criterion is still missing. Lastly, we have 
constructed examples for all the classes and in particular 
for tripartite entangled state with positive partial trans- 
pose. Using the separability criterion for multi-mode bi- 
partite Gaussian states |23] the results presented above 
can be extended to cover the case of n modes at location 
C. Nothing changes in the argumentation to distinguish 
3-party biseparable from fully separable states [the addi- 
tional modes are taken care of automatically in Eqs. (^^] . 
However, the separability criterion of [p2[ is now neces- 
sary to determine the properties under bipartite splitting, 
since for AB-C we deal with a 2 x n state and PPT is 
then no longer sufficient for biseparability pCl]. 



It is worth pointing out that the separability crite- 
rion can be checked experimentally. The CM 7 can be 
measured, and thus the criterion is entirely formulated 
in terms of quantities that are measurable with current 
technology. 

Gaussian CV states promise to be a fruitful testing 
ground for quantum nonlocality: Pure entanglement is 
comparatively easy to create in quantum optical exper- 
iments, as described in ||2^ . Likewise, tripartite bound 
entangled states are experimentall y acce ssibl e: the states 
discussed in the examples Subsec. [VB and IV C can be 



obtained by mixing differently displaced pure Gaussian 
states. 

The study of entanglement of multi-party Gaussian 
states is still in a very early stage. For example, no 
work has to our knowledge been done on the interest- 
ing cases of more parties and modes. But even for the 
simple three-mode case there are important open ques- 
tion. In particular nothing is known about the distilla- 
bility of tripartite states. As in Rcf. for qubits, it is 
easy to see that Gaussian states in Classes 3 and 4 can- 
not be distilled at all and are therefore bound entangled. 
For this, we consider N copies of a Class 3 state p, and 
apply an arbitrary local quantum operation Viocc con- 
sisting of a classically correlated sequence of operations 
of the form V = Va ® Vb ® Vc- Since p is in Class 3 
we can write p®^ as a mixture of AB — C product states 
Sfc PkPAB k ® Pc'k ^"^^ ^ mixture of AC — B product 
states J2kPkP^AC k ® P^B^k- After applying an operation 
such as V the resulting state p = V{p®'^) will still be sep- 
arable along these cuts, and no sequence of operations V 
can change this. Thus p is bound entangled. 

Whether all states in Class 2 may be distilled to max- 
imally entangled states between the two non-separable 
parties is an open question. If this were shown, it would 
follow that all states in Class 1 could be distilled into 
arbitrary tripartite entangled states. 
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APPENDIX A: POINTS OF INTERSECTION 



(A6) 



As shown in Theorem |^ a state is separable iff solu- 
tions to Ineqs. ( p^ ) are found among the points of inter- 
section of the curves described by the equalities (p8|), or 
the centers of the three sets. Here we give the formulas 
to directly calculate these points from 7. 

The centers of circle and the ellipses have already been 
shown to be 



(o,or, 

det TV + 1 

fci 

det TV -I- 1 ~ 



(Al) 



where N, N were defined in (|9|) , L in ( P6| ) , and ki , ki 
after (p9|). The intersections of the borders of C,£,£ 
are calculated as follows. Consider first the t wo el l ipses , 
whose borders are defined by the equalities (28b, 28c). 
Dividing by trA^, respectively by tvN and subtracting 
the two equalities we find that a point on both d£ and 
d£ must lie on the straight line Gee defined by 

(det N + 1 + L^C) /tr A = (dct A + 1 + L^^ /trTV, 

(A2) 

where ^ = {y,z). Q^e can be parameterized with s G R 
as gee + sfee, whcrc 



gee = 



det A 1 det A + 1 



trA 



trA 



L'/WL' 



(A3) 



i/trA Q and /eg is a vector or- 



where L' = L/trA 
thogonal to L' . 

Inserting Qee hi the equation ( 281: ) for d£ we obtain a 
quadratic polynomial in s, whose roots s^g (if they are 
real) give the intersection points. For the intersections of 
dC with the ellipses we proceed similarly. In summary, 
we get for the intersection points 



ge 



s/ee ; 



(A4) 



(A5) 



where the vectors g^, x = ce, ce are 

gee = (trAV'Tf+T-detA- 1) (A7a) 

fee is a vector orthogonal to L, and re is the smaller of 
the two radii 



rr = mm 



v/(trA)V4-l, ^(trA)74-l| . (A8) 



<?cej fee are defined likewise for tilded quantities. And, fi- 
nally, by s^g, we denote the real roots of the quadratic 
polynomials 

Pee(s) = {L'^{geS + sfee)+detN+lf - 

{tvNf {l + \\ge, + sfe,\\^) , (A9a) 



- \\gx + sfx\\^,x = ce,ce. 



(A9b) 



Thus all the nine candidates are given in terms of A, A 
which can be directly obtained from 7. 



APPENDIX B: CHARACTERIZATION OF K 



Here we show that ^"(7) as defined in Eq. ( p3^ ) coin- 
cides with the (real) span of the vectors belonging to the 
kernels of 7 + Jxl^^Jx- This fact automatically follows 
from the following 

Lemma 7 (Characterization of K^j)) 

Let f — fn. + if I, where fn, and fi are real. Then 
f e ker(7 - iJ^) iff fj = l^^Jxfn and both fn and fj 
belong to the kernel of j + JxJ~^Jx- 

Proof: Taking the real and imaginary part of the 
equation {'y~iJx)f — we find jfn + Jxfi = and 7// — 
Jxfn = 0. Since 7 must be invertible we obtain from the 
second equation that // = 7"^ Jr/fl. Using now the first 
equation we find that {'~f+JxJ~^Jx)fR — 0. Analogously, 
{1 + Jx"f~^ Jx)fi — 0. The same argumentation holds for 
the other direction of the proof. ■ 
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